Quasi-Spherical Light Cones of the Kerr Geometry 



Frans Pretorius and Werner Israel 



Canadian Institute for Advanced Research Cosmology Program 
Department of Physics and Astronomy 
University of Victoria 
P.O. Box 3055 STN CSC 
Victoria B.C., Canada V8W 3P6 



Abstract 

Quasi-spherical light cones are lightlike hypersurfaces of the Kerr geometry that are asymptotic to 
Minkowski light cones at infinity. We develop the equations of these surfaces and examine their properties. 
In particular, we show that they are free of caustics for all positive values of the Kerr radial coordinate 
r. Useful applications include the propagation of high-frequency waves, the definition of Kruskal-like 
coordinates for a spinning black hole and the characteristic initial- value problem. 



PACS numbers: 0240, 0420 



1 Introduction 



The Kerr geometry is a vacuum spacetime with a Weyl tensor of Petrov type D. According to the Goldberg- 
Sachs theorem |lj], it therefore possesses two congruences (ingoing and outgoing) of shear- free lightlike 
geodesies. Historically, these congruences played an essential role in the discovery of the Kerr solution 
[E|, because the metric takes a simple explicit form in coordinates adapted to them. 

In the original |^] , || (Eddington-Kerr) form of the metric, the ingoing congruence consists of parametric 
curves of the Eddington- Kerr coordinate t>, often referred to as an "advanced time". This nomenclature can 
be misleading, because hypersurfaces of constant v are actually spacelike, not lightlike, if the Kerr rotational 
parameter a is nonzero. (The lightlike congruences are "twisting", i.e., not orthogonal to any hypersurfaces, 
and hence not tangent to lightlike hypersurfaces if a ^ 0). 

Thus, the Eddington-Kerr form of the metric, and the Boyer-Lindquist form Q derived from it by an 
elementary coordinate transformation, correspond to a "threading" [|| (i.e., a one-dimensional foliation) of 
the manifold by twisting lightlike geodesies. However, a "slicing" (three dimensional foliation), in particular a 
light-like slicing, is for many purposes more advantageous and often corresponds more closely to the physics. 

Light cones of the Kerr geometry have not previously been a subject of systematic study to our knowledge. 
Perhaps this stems in part from a feeling that such hypersurfaces would quickly develop caustics because of 
the twist inherent in the metric. Yet as characteristics these surfaces obviously play a key role in the physics 
of Kerr black holes, e.g in the propagation of wave fronts and high-frequency waves, and in characteristic 
initial value problems. 

In this paper we first derive the general solution for axisymmetric lightlike hypersurfaces of the Kerr 
geometry (Sec. |^). We then focus on a particular foliation (invariant under time-displacement) by (either 
ingoing or outgoing) "quasi-spherical" light cones, defined as degenerate 3-spaces whose spatial sections 
are asymptotically spherical at infinity. These hypersurfaces become the standard Minkowski light cones 
when the Kerr mass parameter m vanishes; if a = 0, m =/= 0, they are the familiar spherical light cones, 
t ± r* = constant, of Schwarzschild spacetime (where r» is the "tortoise" coordinate). Quite generally, for 
arbitrary m > and a, they have the remarkable property of being free of caustics for all positive values of 
the Kerr radial coordinate r. 

These appealing features suggest that coordinates adapted to quasi- spherical surfaces should be especially 
convenient and useful for the description of the Kerr geometry. The drawback (and it is considerable) is that 
the metric is no longer expressible in an explicit elementary form, because the new coordinates are elliptic 
functions of the Boyer-Lindquist coordinates. 

We now briefly outline the contents of the paper. We begin in Sec. 2 by solving the eikonal equation 
for axisymmetric lightlike hypersufaces of the Kerr geometry. The general solution is expressed in terms of 
integrals of elliptic type, and involves an arbitrary function / of one variable, to be fixed when the "initial" 
(e.g. asymptotic) shape of the surface is given. (If m — the geometry becomes flat and the integrals 
can be reduced to elementary form; in Sec. 3 we digress briefly to illustrate this.) Some general properties 
of axisymmetric lightlike surfaces are developed in Sec. 4. In Sec. 5, guided by the results of Sec. 3 for 
Minkowski light cones, we make the special choice of / appropriate for asymptotically spherical light cones in 
Kerr space. The key result, that these are free of caustics for r > 0, is demonstrated in Sec. 6. In Sees. 7 and 
8 we show how the "quasi-spherical" coordinates r* , 6** defined by these light cones can be used to construct 
high-frequency solutions to the wave equation and Kruskal-like coordinates for Kerr black holes. Sec. 9 
summarizes the limiting cases that can provide serviceable approximations for the elliptic functions which 
appear in the quasi-spherical form of the metric. Finally, in Sec. 10 we present some results of numerical 
integrations for the inward evolution of the light cone generators. 

2 Axisymmetric lightlike hypersurfaces 

The equation of an arbitrary axially symmetric lightlike hypersurface of the Kerr geometry is expressible in 
terms of elliptic integrals as we now proceed to show. 

We write the Kerr metric in its standard (Boyer-Lindquist) form 
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,o S , 9 „ ,„o „o 9„,t 4mar sin 2 6 , , , ( 2mr\ , 
ds 2 = —dr 2 + Y,d9 2 + R 2 sin 2 9d(/> 2 d(f>dt - ( 1 - — J * 2 , (1) 

where 

S = r 2 +a 2 cos 2 0, i? 2 = r 2 + a 2 + 2ma ^ sin2 g | A = r 2 + a 2 - 2mr, (2) 
and we note the useful identities 

Si? 2 = (r 2 + a 2 ) 2 - Aa 2 sin 2 9, (3) 
5005tt - 9% = ~ A sin2 S - ( 4 ) 

The equation 

v(t, r, 8) = t + er*(r, 9) = constant, e = ±1 (5) 

represents an axisymmetric lightlike hypcrsurface (ingoing for e = 1, outgoing for e = — 1) if g a @ {d a v)(dpv) 
0, i.e., if 

A(d r r*) 2 + {der*) 2 = {r 2 + a 2 ) 2 /A - a 2 sin 2 6. (6) 

It is easy to obtain a particular (separable) solution of (||) depending on two arbitrary constants (a 
"complete integral") by adding and subtracting an arbitrary separation constant a 2 A on the right -hand 
side. Let us define 

P 2 (M) = a 2 (A-sin 2 6>), Q 2 (r, X, m) = (r 2 + a 2 ) 2 - a 2 XA. (7) 
(Note the useful identity 

Q 2 + AP 2 = Si? 2 , (8) 
which follows from (^).) A complete integral 

r* = p(r, 6, A, m), (9) 
of (0) is then obtained by integrating the exact differential 

dr* = (Q/A)dr + Pd6 (10) 



at fixed A. When ( |10[) is integrated, a second, additive integration constant appears which we shall denote 
as et 2 /(A)/2, where / is an arbitrary function. 

We next proceed in the usual way to promote this complete integral, depending upon the arbitrary 
constants A and /(A), to a general solution involving an arbitrary function. 

In ([)]), p is a function of three independent variables r, 9 and A (not counting to and a), and a more 
complete expression for its differential is 

dp=(Q/A)dr + Pd9+(a 2 /2)Fd\, (11) 

where F is the partial derivative 

(a 2 /2)F(r, 9, A, to) = 9 A p(r, 9, A, to). (12) 
Its explicit form may be taken to be 
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r d9' r°° dr' 

Up to now, we have taken A to be a constant, i.e., dX = in (|ll|). But we achieve the same effect (i.e., 
il]) still reduces to (|Io| ) ) even when dX ^ provided we require that F = 0. In other words: the function 
t (r, 9), given by @ with A now a function A(r, 0), remains a solution of (& provided its extra dependence 
on r, 9 through A does not change the algebraic form of the differential ( |10|) . This will indeed be so if we 
impose the constraint 

F(r,6,\,m) = 0. (14) 

This condition fixes the dependence of A on r, 9 for any given choice of /(A). Thus we now have a general 
solution, depending upon an arbitrary function f[X(r,9)]. 
The explicit form of the general solution M) is then 



+ / d0'P(9',X) + la 2 f(X), 



where the radial dependence has been arranged to ensure convergence of the definite integral at its upper 
limit. When performing the integrations in (|l^) and Jl5|), we are allowed to treat A as merely a passive 
constant parameter - it is in fact a function of the limits of integration r,9, not of the integration variables 
- with the constraint ( p"4| ) imposed a posteriori. It is possible, though not especially illuminating or useful, 
to express each of these integrals in terms of standard elliptic integrals. 

3 The case m = 0: light cones in Minkowski space 

As a simple illustrative example (and because we shall need some of the results later), we specialize in this 
section to the case m = 0. The Kerr line-element reduces to the metric of flat spacetime expressed in terms 
of oblate spheroidal coordinates r, 8 - related to Cartesians by 

x 2 + y 2 = (r 2 + a 2 ) sin 2 9, z = rcos9. (16) 

The second integral in (Q) can be reduced to the same general form as the first integral when m = 0. 
Assume that we are in a domain where A(r, 9) < 1. Setting 



Qo(r', A) = Q(r', A, m = 0) = Vr' 2 + a 2 ^r' 2 + a 2 (l - A), (17) 
and making the substution 



Vr^Asin X 7V A - sin X', (18) 



we find 



dr'/Q = dx'/a^X- sin 2 x > . (19) 
Thus the two integrals in ( ^3|) can be combined into a single integral, to give 



„F[r.9.\.,n = ())= I dX [ +g'(A), (20) 

*(r,A) y/X - sin x' 
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where 

g(X) = af(X) + 2 f ' Ja- sin 2 (21) 
Jo 

and the form of x( r > A) is given by the "unprimed" version of (fj"q); equivalently, 



tanft^ = tanx, A = sin 0*. (22) 

r 

The function g(X) is arbitrary. As an example, let us consider the simplest choice, g(X) = 0. The 
functional dependence A(r, 9) corresponding to this choice is determined by the constraint F = 0, which 
requires that x{ r -> X) = 9 by inspection of (|20|). Thus (E2) gives 



tan ft, = ^!±Z tan = (23) 

r z 

by (|l6|), showing that ft* is the spherical polar angle. With A(r, 0) known from (|3|), it is straightforward 
to integrate (nffl to obtain 



■''* 



\Jr 2 + a 2 sin 2 ft (24) 

Thus r* is the usual spherical radius, and our solutions v = t ± r* = constant of (^) are in this case the 
Minkowski light cones. 

4 Axisymmetric lightlike hyper surfaces: general properties 

Returning to the general case, we shall now derive a number of results applicable to all axisymmetric lightlikc 
hypersurfaces of Kerr. 

Since the function A(r, 9) is determined by the constraint F = 0, the partial derivatives of A can be read 
off from dF = (pt): 



lidX = -dr/Q + d6/P, n = -dF/dX. (25) 
It follows from ( |l0| ) and ( p5| ) that Vr„ • VA = 0, i.e. that r* and A are orthogonal with respect to the 



intrinsic 2-metric 



da 2 = Y,(dr 2 /A + d9 2 ) (26) 



of the spatial sections ((f), t) = constant of the Kerr geometry. Since A is independent of <f> and t, it follows 
further that A is constant along the lightlike generators, i.e., 

l a d a X = 0, i a = -d a v = + en). (27) 
If r* and A are adopted as coordinates in place of r and 0, the 2-mctric ( p6|) becomes 

do- 2 = R- 2 (Adr* 2 + L 2 dX 2 ), L = fiPQ, (28) 
where we have made use of @, © and (|). The identity (@) can be re-arranged in the form 

l-^)+-|i? 2 sin 2 0= A, (29) 



where 



qm 2mar 

c B = -^ = __ (30) 
g^ hn z 
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is the Bardeen or ZAMO angular velocity d<p/ dt which characterizes orbits having zero angular momen- 
tum. This enables us to recast the Kerr metric (Q) in the form 



Thus 



A T 2 

ds 2 = —,{dr 2 - dt 2 ) + —d\ 2 + R 2 sin 2 9{U - uj B dt) 2 . (31) 
R z Rr 



g r * r * =R 2 /A = -g tt , (32) 



from which the null character of the 3-spaces t ± r* = constant is directly evident. The degenerate 
intrinsic metric of these 3-spaces is 

{ds 2 ) LL = (L/R) 2 d\ 2 + i? 2 sin 2 6>(d</> - uj B dt) 2 , (33) 

showing that the generators rotate with ZAMO angular velocity relative to stationary observers at 
infinity. (This is directly obvious from (p7p, which shows that 1$ — for an axisymmetric hypersurface). 
From (^) we see that caustics will develop when the degenerate volume element tends to zero, i.e. when 
L sin 9 — > (recalling that A has a fixed value along each generator) . 

The integrability of (^5|) for the exact differential d\ provides a condition on the integrating factor u. 
This takes the form of an evolution equation for u along the generators. A straightforward calculation, which 
calls upon (g), (@), © and @ itself, yields 

(W*h = Tl£&> (34) 

where the subscript A indicates that the partial derivative is being taken at fixed A. 
Inversion of the differential relations (|l^) and (25) gives 

ZR 2 dr = AQ(dr*- uP 2 d\), ZR 2 d6 = P(Arfr„ + uQ 2 d\), (35) 
which will be of use in subsequent sections. 



5 Quasi-spherical light cones 

From here on we confine attention to a special class, "quasi-spherical light cones" , defined as lightlike 
hypersurfaces which asymptotically approach Minkowski spherical light cones at infinity. Define an angle 
0*M)by 

A = sin 2 9*. (36) 

In Minkowski space, the generators A = constant of light cones are radial straight lines, suggesting (as 
we confirmed in Sec. ||) that 9* is the spherical polar angle. For r >> a, the oblate coordinates (r, 9) become 
indistinguishable from spherical coordinates, so we have the condition 

9*{r = <x>,6)=6 (37) 

for spherical light cones in Minkowski space. This asymptotic condition must therefore also hold for 
quasi-spherical surfaces in Kerr space. This fixes the arbitrary function /'(A) in (|l3|). We conclude that 

f°° dr' f e * d& 

F(r,6,\m)= - f\ - / (38) 
J r Q{r',X,m) J e P{9',\) 

generates the solution for quasi-spherical surfaces. The radial function p(r, 9, A, m) is given by (|l^) with 
(compare (||) with 5(A) = 0) 
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a/(A) = -2 / * y/X-sin 2 6'd6'. (39) 
Jo 

The equation of the hypersurface is then v = t ± (r, 0) = constant, with r* = p(r,9, X(r,9)) and the 
function A(r, 6*) = sin 2 0* determined by the constraint F — 0. 



6 No caustics for positive r 

We now prove the rather remarkable result that quasi-spherical light cones are free of caustics for all positive 
values of the Kerr radial coordinate r. 

This is trivially true if m = 0, when these surfaces are simply light cones in Minkowski space with vertices 
at the spatial origin, represented by r = 0, 9 = or n in oblate spheroidal coordinates according to (|l6|). We 
shall prove that it is true a fortiori for to > by effectively showing that when to is larger, the generators 
converge less rapidly as r — > + . 

As noted in Sec. ^, formation of a caustic along a generator is signalled by 

fj,PQsm6^>0. (40) 

We consider in turn the behaviour of each factor in ( fiof ) along an ingoing generator, so that r is pos- 
itive and decreasing, with A fixed. Because of the equatorial symmetry we need only consider "northern" 
generators, i.e., we may assume (which specifies the initial, asymptotic value of 9 when r = +oo) to be 
acute, and P(6, A) positive, at least initially. (Note from ( |To|) that equatorial symmetry of r*(r, 6) implies 
P(A,tt-0) =-P(A,0).) 

Since 9 decreases with r at fixed A according to (p5|), the factor P = a — sin 9 must increase and 
remain positive. Any possible caustic in the Kerr positive-r sheet cannot arise from the behaviour of P. 
We proceed to consider the other factors in (ff"o|). From the definition (0), 

Q(r, A, to) > Q{r,X,m = 0) > (r>0,TO>0). (41) 
We shall further show that 



(d9/dm) r:X > and (dfj,/dm) r , x > (r > OA < tt/2). (42) 

Thus, none of the three factors u, Q and sin# can reach zero sooner for positive to than they do in flat 
space, and hence they do not reach zero for any positive value of r. 

To establish (ff2|), we note that the function 9(r, A, to) is determined by the vanishing of F(r, 9, A, to) as 
given by (38). Taking the differential at fixed r,A and using (R) gives 

poo 

{d9/dm) r:X = a 2 \PF I(r,X,m)= r'dr' /Q 3 (r' , A, to), (43) 

J r 

which is manifestly positive. Turning to /i, this is defined as a function of r,9,X,m by [i = —dF/dX 
according to (B5|) . Inverting the order of partial differentiation and using 



du d dF d 2 N r,A\ 
= ^( flA/ )' ( 44 ) 



dm dX dm dX 
du d dF d 1 a 2 
~dl = ~£)X~d9 = ~dXP = 2P^ 



(45) 



H6nc6 

du\ dn du ( d9\ of la 2 A\ r 3 4l T 



G 



where 



J(r,X,m)= r'A(r')dr'/Q 5 (r',\,m). (47) 



Each term in (Q) is manifestly positive. This completes the proof. 

The fact that adding mass to the Kerr field should reduce the focusing of ingoing radial light rays is 
at first sight paradoxical, but can be explained by the circumstance that the material source is not at the 
origin, but located on the singular equatorial ring x 2 + y 2 = a 2 , z = (see (|X6|) ) . The source has a peculiar, 
"demi-pole" structure [||, made possible by the double-sheeted structure of the Kerr manifold: a ring of 
positive mass in the sheet r > is bonded to a ring of negative mass in the sheet r < 0. Light rays heading 
inwards toward the origin in the sheet r > are deflected outward by the ring and defocused. Repulsive 
effects have become dominant by the time they pass through the disk r — : the rays are then refocused 
and form a caustic in the sheet r < 0. This description is more than just hand-waving: the Keres Newtonian 
analogue model of the Kerr field [^| has a source structure and equations of motion closely resembling Kerr 
(apart from frame dragging), and displays precisely this behaviour. 

Figure [l] (a result of numerical integrations described in Sec. 10) shows the ingoing generators mapped 
onto the flat background of the Kerr-Schild decomposition, using the rectangular coordinates defined in (|l^) . 
(Effects of frame dragging are not shown). According to the Keres-Kerr model, the gravitational "force" 
should become repulsive for r ~ a, and, indeed, the generators have inflection points near this radius. 

7 Eikonal solution of the wave equation 

Since the quasi-spherical light cones £±r* = constant are characteristics of the wave operator, the coordinates 
r„, A are well adapted for representing asymptotically spherical high-frequency modes. 
The wave equation for \I/(r„, A, (j>, t) on the Kerr background ( |3l| ) is 

= -d 2 )^ + {drjni){dr,^) + i- 1 dx{ 1 - 1 l\siii 2 ed^)+ (48) 

where 7 = [iPQsmO gives the area element (degenerate volume element) on the light cone (|3 
Introducing the ansatz 



* = A)e i " (t+er * ) e^ (e = ±1) (49) 



into (E8h, we find 



where 



2eiu; 



<9 r „$ + $ Q^r, In 7 - einujBj 



B = 0, (50) 



B = dl$+ {d r , ln 7 )(0 r .$) + 7" 1 9a (7" 'A sin 2 9d x <S>) - „ 9 . 2 „ 1 - — U (51) 

R^sm \ 2-> J 

has no explicit w-dependence. We can re-express in ( |50"| ) as a partial derivative with respect to r*. 
Defining an angular function 

f r 2mar' 

aM) =/ wm^) dr (52) 



and recalling (p5|),we find 
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Figure 1: Null generators of the quasi-spherical light cones projected onto the Cartesian plane for m = 1, 
a = 1/2. The azimuthal twist (not shown in this two-dimensional representation) depends on the choice of 
azimuthal coordinate; it is zero for the coordinate (j) e introduced in Sec. 8. In order to follow the generators 
through the equatorial disk r = and into the negative-r sheet, the "southern half" (9 > n/2) of the 
positive-r sheet has been cut away, and replaced by the "northern" half (9 < ir/2) of the negative-r sheet. 
Note that only well within the outer horizon does the mass of the spacetime significantly affect the r, 9 
behaviour of the generators (which would continue as straight lines intersecting the origin if m were zero). 
The caustic emanates from the singularity r = 0, 9 = tt/2 (if m > 0) and creeps inwards toward the axis 
and "downwards" toward increasingly negative values of r. 
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da \ 2mar ( dr \ 2mar 

uj B - (53) 



dn) x AQ \dnj x Si? 2 
Therefore (|50D can be rewritten 



2eiLud rt (VJ e ~ ema ®) = -^e- tlna B. (54) 

If lo is large, we can apply an iterative procedure to ( |54"|) and (|5l]) to develop the solution for <E> in inverse 
powers of iv. In the lowest approximation, we simply equate the coefficient of uj in (|54| ) to zero. This yields 
the "eikonal approximation" 

» 7 - 1 /2y( A ) e ^(t+er») e in(<H+e Q! ) ^ 

with an arbitrary function /(A). By linear superposition we can build from this an arbitrary high- 
frequency solution 

^ 7 -i/2F{\,(j> + ea,t + en), (56) 
where F is an arbitrary function of its three arguments varying rapidly with time. 

From ([55]) and ( |53] ) we see that crests of high-frequency azimuthal ((^-dependent) waves twist about the 
axis with angular frequency ujb, the ZAMO angular velocity. Because lob and a depend on A, propagation 
of azimuthal waves produces latitude-dependent phase shifts. The latitude-dependence is, however, always 
bounded, even when a — > ±00 at horizons. 

In particular, wave-tails propagating inwards from the event to the Cauchy horizon of a spinning black 
hole experience a blucshift which is modulated by an oscillatory, latitude-dependent factor when n =/= - 
an effect first noted by Ori [Q. However, since \n\ is no larger than the multipole order £, this effect is 
more than offset by the natural power- law decay of the higher multipole wave-modes with advanced time 
v, \& ~ v -( 2e + 2 ) . It does not affect the uniformity of the leading- order divergence e^ K ~^ v of blueshift at 
the Cauchy horizon, where k_ (the inner-horizon surface gravity) is a constant. These features, which are 
expected to extend at least qualitatively to waves having (initially) lower frequencies, are important when 
considering the back-reaction of blueshifted wave-tails on the geometry near the Cauchy horizon §| . 



8 Kruskal coordinates 

It is straightforward to transform the Kerr metric ( |3i"| ) into a Kruskal-like form. We introduce retarded and 
advanced times u and v, and associated Kruskal coordinates U and V, by the definitions 

dU , , . dV , ,. . . 

— = du — dit — r*J, — — = dv = d(t + r* ) . (57) 

ku kv 

Here, k is the surface gravity of the horizon under consideration, defined for the outer (r — r + ) and inner 
(r = r_) horizons by 



K± = ±\/in 2 -a 2 /(r±+a 2 ). (58) 

Then ([n]) becomes 

1 A L 2 

ds 2 = —— — dUdV + —dX 2 + R 2 sin 2 9(d0 - uj B dt) 2 (59) 
k z R z U V H z 

with 

-UV = e 2Kr * and - V/U = e 2Kt . (60) 
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The first (dUdV) term is manifestly regular at the horizon sheets U — and V = 0. The last term, 
involving the Boyer-Lindquist coordinates <j> and t, is not. We therefore define the advanced and retarded 
angular coordinates <f>+, 4>- by 

<j> e =(j) + ea{r,X) (e = ±1) (61) 
where a was defined in (p3). It is straightforward to show that 

da = uj B dr* - Nd\, N = ^P 2 uj b + ma s I (62) 
with I defined as in (|43|). Thus the last term in ( |59| ) can be expressed in a manifestly regular way as 

d<f> — LUsdt = d(p e — u)Bdu e + eNdX, (63) 

with u+ = v, U- = u. In ( pq ) we can choose either sign for e, depending on which sheet of the horizon 
is of interest. For example, </>+ and it+ are regular on the "upward" (future) sheets of both outer and inner 
horizons, with the nice bonus feature that u+ is constant over each ingoing light cone and </>+ is constant 
along each ingoing generator. 

Is there a single azimuthal coordinate which regularizes the metric simultaneously on both past and 
future sheets of (say) the outer horizon, including the bifurcation surface? Since ojb takes a constant value 
u>H = a/2mr + over this horizon, a possible choice for the desired coordinate is $ = — Wffi, and <!> in 
fact agrees with <fi + on the future sheet and 4>- on the past sheet. But $ develops the undesirable features 
characteristic of rigidly rotating axes in the outer regions of the space. 

We briefly mention that following standard methods one could construct a Penrose diagram of the Kerr 
spacetime using the Kruskal-like coordinates (|6^), valid up to formation of the caustic surface. Such a 
diagram would not look any different from the usual textbook examples || except that a single 2-dimensional 
diagram is all that would be needed to illustrate the causal structure of the spacetime. This is because the 
ingoing (V — constant) and outgoing (U = constant) lightlike congruences intersect at surfaces of constant 
r», represented by a single point on a diagram where the compactified coordinate system is derived from 
(pOh . The causal future of observers on the surface of intersection is entirely contained within the future- 



directed wedge of U = constant, V = constant (as is evident from (31) by noting that when |<ir*| > \dt\ the 
line-element is spacelike). 



9 Limiting cases and approximations 

The quasi-spherical coordinates r* and A = sin 2 9* , and the coefficients in the quasi-spherical forms ( |3l| ) and 
([39]) of the Kerr metric are complicated elliptic functions of the Kerr-Boyer-Lindquist coordinates r,9. It 
may therefore be useful to record here the simple forms these functions reduce to in the two opposite limiting 
cases, m/a — > and m/a — > 00. 

(i) m = 0, a =/= 0: This was studied in Sec. 3, and we simply list the results. 



r* = V r 2 + a 2 sin 2 6, tanfl* = r + " tan g ; (54) 

r 

r(r 2 + a 2 ) ^ a 2 sin#cos# „ 1 sin cos ,„_. 

P = r=„ „ =^T> » P = o 751 a — a ,2 > ( 65 ) 



Vr 2 + a 2 sinV r 2 + a 2 sin 2 6»' 2 ( sin °* cos 9 

1 r 2 + a 2 sin 2 6 . . 

7=x 7, ■ (66) 

' 2 cos0* x ' 

(ii) m/0, a = : This is just Schwarzschild spacetime, and one readily finds 



/t 



dr 1 
— 2m/ r sm20 
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' 2cos0' v ' 

Comparison with numerical results described in the next section suggest that the m = 0, a / case 
functions can provide serviceable analytic approximations outside of the outer horizon for broad ranges of 
m and a. For example, with m = 1 and a = 1/2, all of the quantities except r* in (|64|), ( |65| ) and ( |66| ) differ 
by no more than roughly 1 to 3 percent from the true (numerically integrated) values, outside of r+. (r* 
converges only logarithmically to the Minkowski value in the limit r — > oo, and of course diverges at the 
horizon) . 



10 Numerical evolution of the generators 

In Figure [TJ we illustrated the behaviour of the generators and the nature of the caustic that forms in the 
negative r sheet of the Kerr manifold. The A = constant curves were obtained by numerically integrating 
the evolution equations of the null generators of the hypersurface: 

£ a = dx a /dr = -goPdpv, (69) 
with affine parameter r. In particular 

f = dr/dT = -eQ/Y,, 9 = dQ/dr = —eP/Y,. (70) 

To detect the formation of a caustic (see ^) we track the evolution of /i, and to reconstruct surfaces 
of constant r* in the r, 9 plane we calculate the change of r* along the generators. If we treat r as the 
independent variable, from (fro|), ([Tsl), (|25| ) and ( |34|) we find 

(dd/dr) x = ~, (71) 

{dn/dr) x = Si? 2 /AQ, (72) 

and 

= (73) 

At the horizons r* and its derivative with respect to r diverge. These divergences are coordinate singu- 
larities which can be avoided numerically by subtracting off the infinities that occur there. Thus define f 
as 



ln\r — r + \ ln\r — r_| 
2k_i_ 2k- 



(74) 



where is the surface gravity at the outer(r + )[inner(r_)] horizon as defined in (f58[). So we actually 

integrate f along the generators and retrieve r* from the result using (Q). 

For asymptotic intitial conditions for a and 9, we use the relations listed in ( |64| ) and ( |65| ) for a sufficiently 
large initial r. Using p^ ) and ( ^6| ) one can show that (i(r, \,m) — u(r, \,m = 0) « (a 2 Am/8sin#cos0)(l/r 2 ) 
and 9(r,\,m) — 9(r,X,m — 0) ~ (a 4 Amsin0cos0/4)(l/r 5 ) for large r along a given generator. The initial 
condition for r* is arbitrary (r* + constant is still a solution to @). 

Figure ^ below is a close-up within the inner horizon of the case illustrated in Figure [jj showing the 
projections of surfaces of constant r* onto the -y/a; 2 + y 2 , z plane. This is the only region where r* = constant 
surfaces start showing significant deviations from the spheres of the massless scenario. Note that the caustic 
surface does not coincide with a hypersurface slice r* = constant: following the generators inward a caustic 
first develops where the hypersurface meets the ring singularity, after which it quickly "unravels" . We only 
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Figure 2: Null generators and surfaces of constant r* projected onto the Cartesian plane (m — l,a — 1/2). 
This is a close-up of Figure [l] focusing on the region between the inner horizon and caustic. As in Figure 1 
the northern half of the negative-r sheet replaces the southern half of the positive-r sheet to illustrate the 
evolution of the generators through the equatorial disk. 

present plots for m = 1 and a — 1/2; except for variations in scale there arc no qualitative differences in the 
shapes of the curves and surfaces for arbitrary non-zero a and positive m. 

11 Concluding Remarks 

Despite the complexity of the coordinate transformations linking them to the familiar Boyer-Lindquist coor- 
dinates, the quasi-spherical coordinates r*, 8* introduced in this paper and the light cones associated with 
them provide new and useful insights into the structure of the Kerr geometry and wave propagation in Kerr 
spacetime. We anticipate that they will find increasing use as their special advantages become apparent. 
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